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Abstract. In this paper we prove a new matrix Li-Yau-Hamilton estimate for Kahler-Ricci 
flow. The form of this new Li-Yau-Hamilton estimate is obtained by the interpolation con- 
sideration originated in [Chi]. This new inequality is shown to be connected with Perelman's 
entropy formula through a family of differential equalities. In the rest of the paper, We show 
several applications of this new estimate and its linear version proved earlier in [CN]. These 
include a sharp heat kernel comparison theorem, generalizing the earlier result of Li and 
Tian, a manifold version of Stoll's theorem on the characterization of 'algebraic divisors', and 
a localized monotonicity formula for analytic subvarieties. 

Motivated by the connection between the heat kernel estimate and the reduced volume 
monotonicity of Perelman, we prove a sharp lower bound heat kernel estimate for the time- 
dependent heat equation, which is, in a certain sense, dual to Perelman's monotonicity of 
the 'reduced volume'. As an application of this new monotonicity formula, we show that the 
blow-down limit of a certain type III immortal solution is a gradient expanding soliton. In 
the last section we also illustrate the connection between the new Li-Yau-Hamilton estimate 
and the earlier Hessian comparison theorem on the 'reduced distance', proved in [FIN]. 



§0 Introduction. 

In [LY], Peter Li and S.-T. Yau developed the fundamental gradient estimates for positive 
solution u(x,t) to the heat equation (J^ — A)u(x,t) = 0. On a complete Riemannian 
manifold M with nonnegative Ricci curvature, they also derived a sharp form of the classical 
Harnack inequality (cf. [Mo]) out of their gradient estimates. Later in [H2], Richard 
Hamilton extended the estimate of Li- Yau to the full matrix version on the Hession of 
logu, under the stronger assumption that M is Ricci parallel and of nonnegative sectional 
curvature. More recently, in [CN], H.-D. Cao and the author observed that if M is a 
Kahler manifold with nonnegative bisectional curvature, one can obtain the Hamilton's 
matrix version estimate on the complex Hessian of log u without the assumption of Ricci 
being parallel. Following [NT1], We called our estimate in [CN] a Li-Yau-Hamilton estimate 
(or inequality, which is also referred as differential Harnack inequality in some literatures) . 
For Ricci flow (Kahler-Ricci) flow there also exists the fundamental work of Hamilton 
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[HI] (H.-D. Cao [CI]) on the matrix Li-Yau-Hamilton inequality for curvature tensors. 
See also [CC1-2], [CH], [CK1], [NT1], etc, for the later developments, [A] and references 
therein for the Li-Yau-Hamilton estimates of other geometric flows. The relation with the 
monotonicity formulae was discussed in [N2, N4]. 

In this paper we shall prove a nonlinear version of the Li-Yau-Hamilton estimate of 
[CN], for time dependent Kahler metrics evolving by Kahler-Ricci flow. The new matrix 
inequality asserts that if (M, g(t)) is a complete solution to Kahler-Ricci flow -§j:g a p(x, t) = 
—R a p(x,t) with (bounded, in case M is not compact) nonnegative bisectional curvature, 
and if u is a positive solution to the forward conjugate heat equation: ( — A — 71) u(x, t) = 
0, where 1Z is the scalar curvature, then 

u 

(0.1) u a - p + -g a p + uR a - p + u a Vp + u- p V a + uV a Vp > 

for any (1,0) vector field V. The form of this new matrix Li-Yau-Hamilton type estimate is 
a natural one after we found it through an interpolation consideration which was originated 
in [Chi] by Ben Chow. We shall illustrate this interpolation consideration further in the 
following paragraph. 

In [CH], Ben Chow and Richard Hamilton proved a linear trace Li-Yau-Hamilton in- 
equality for the symmetric positive definite 2-tensors evolved by the (time-dependent) Lich- 
nerowicz heat equation, coupled with the Ricci flow, whose complete solution metrics have 
bounded non-negative curvature operator. This result in particular generalizes the trace 
form of Hamilton's fundamental matrix Li-Yau-Hamilton estimate on curvature tensors, for 
solutions to Ricci flow in [HI] . Later, in [Chi] , Chow discovered a very interesting interpo- 
lation phenomenon. Namely he shows a family of Li-Yau-Hamilton estimates in the case of 
Riemann surfaces with positive curvature such that this family connects Li-Yau's estimate 
for the positive solutions to a heat equation with Chow-Hamilton's linear trace estimate on 
the solutions to the Lichnerowicz heat equation, coupled with the Ricci flow, in the case of 
Riemann surfaces. Seeking the analogue of such interpolation in higher dimensions turns 
out to be fruitful. Even though the interpolation itself has not been found directly useful 
in geometric problems, it does play a crucial role in discovering new (useful) estimates. For 
example in [N4] , such consideration led the author to discover a Li-Yau-Hamilton estimate 
for the Hermitian-Einstein flow (what proved there is more general). This estimate of [N4] 
is (crucial) one of the new ingredients in obtaining the (sharp) estimates on the dimension 
of the spaces of holomorphic functions (sections of certain line bundles) of the polynomial 
growth. (For more details, please see [N4], as well as [CFYZ].) The Li-Yau-Hamilton in- 
equality proved in [N4] can be interpolated with the earlier one proved by L.-F. Tarn and 
the author in [NT1]. This interpolation consideration further suggests that there should 
be a one-one correspondence, for the Li-Yau-Hamilton estimates, between the linear case 
and the (nonlinear) case with Ricci flow. Seeking the linear correspondence of the linear 
trace Li-Yau-Hamilton inequality for Kahler-Ricci flow, proved in [NT1], led the correct 
formulation of the Li-Yau-Hamilton's inequality in [N4] for the linear case. On the other 
hand, looking for the nonlinear version of matrix Li-Yau-Hamilton estimate for the linear 
equation in [CN] leads us to formulate the correct form of the (nonlinear version) matrix 
Li-Yau-Hamilton inequality in this paper for the Kahler-Ricci flow. The proof of this re- 
sult is applying the tensor maximum principle of Hamilton, which can not be completed 



LI-YAU-HAMILTON INEQUALITY 



3 



without the generous help from Professor Ben Chow on a certain crucial step. We want to 
record out gratitude to him here. 

As in the most other cases, the new estimate proved in this paper is sharp since it 
holds equality if and only if on expanding Kahler-Ricci solitons. Its proof also makes use 
of the earlier fundamental Li-Yau-Hamilton estimate of H.-D. Cao [CI] (see also [C2]). 
The Riemannian version of (0.1) suggests a new matrix differential inequality on curvature 
tensors, which is different from Hamilton's one in [HI]. Please see section 5 for details. This 
new expression also vanishes identically on expanding solitons. Unfortunately, we have not 
been able to verify this new matrix estimates at this moment. (Please see Remark 5.2 for 
details.) 

A little surprisingly, this new Li-Yau-Hamilton estimate for the Kahler-Ricci flow can be 
shown to be related to Perelman's celebrated entropy formula for the Ricci flow [P], at least 
for Kahler case. Again this is done through Chow's interpolation consideration. Namely, 
one in fact can obtain a family of pre-Li-Yau-Hamilton equalities (a notion suggested to us 
by Tom Ilmanen) such that at one end, the trace of this pre-Li- Yau-Hamilton equality gives 
the Perelman's entropy formula after integration on the manifold, and at the other end 
one obtains the Li-Yau-Hamilton inequality of this paper by applying the tensor maximum 
principle of Hamilton. At the midpoint one can obtain both the entropy formula for the 
linear heat equation proved in [N3] and the Li-Yau's estimate in [LY] for the linear heat 
equation. (The pre-Li-Yau-Hamilton equalities was proved earlier by Ben Chow in [Ch2] 
(see also [CLN]) for the backward Ricci flow on Riemannian manifolds soon after the proof 
of (0.1).) This connection between (0.1) and Perelman's entropy formula suggests that 
the new Li-Yau-Hamilton estimate proved here may be of some fundamental importance. 
Please see Section 3 for the detailed exposition on this interpolation between Perelman's 
entropy formula and the new Li-Yau-Hamilton estimate. One should also refer to the recent 
beautiful survey [Evl] and the stellar notes [Ev2] by Evans for the relation between entropy 
and the Harnack estimates for the linear heat equation (see also [N3] for a different relation). 
[Evl-2] also contain many other applications of entropy consideration in the study of PDE. 

The rest of the paper is on applications of this new estimate, as well as the corresponding 
linear one of [CN]. The immediate consequences include the monotonicity of a quantity 
which is called Nash's entropy, a Perelman type monotonicity (or Huisken type in the 
linear case) of the 'reduced volume of analytic subvarieties' and a sharp form of Harnack 
estimate for positive solutions to the forward conjugate heat equation. It also can be applied 
to proved a sharp heat kernel comparison theorem for any subvariety in complete Kahler 
manifolds with nonnegative bisectional curvature. This generalizes a previous result of 
Peter Li and Gang Tian [LT], in which the authors proved the sharp comparison on heat 
kernels of algebraic manifolds, equipped with the induced Fubini-Study metric (also called 
Bergmann metric in [LT]) from the ambient P m . More precisely, we proved the following 
result. 

Let M be a complete Kahler manifold with nonnegative bisectional curvature. Let 
H(x, y, t) be the fundamental solution of the heat equation. Let V C M be any complex 
subvariety of dimension s. Let Ky(x, y, t) be the fundamental solution of heat equation on 
V. Then 



(0.2) 



K v (x,y,t) < (ivt) m s H(x,y,t), for any x, y G V. 
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The equality implies that V is totally geodesic. 

A more involved application is the localized monotonicity formula and an elliptic 'mono- 
tonicity principle' for the analytic subvarieties. The later leads to a manifold (curved or 
nonlinear) version of Stoll's characterization on 'algebraic divisors'. This localization uses, 
substantially, the beautiful ideas from the the study of mean curvature flow of Ecker in 
[El-2]. We shall give a brief sketch on these results below. 

Let V be a subvariety of complex dimension s as above. Denote by -4y,x (p) the 2s- 
dimensional Hausdorff measure of V D B Xo (p) . Here B Xo (p) is the ball (inside M) of radius 
p centered at x . The elliptic 'monotonicity principle' states the following. 

There exists C = C(m, s) such that for any p' G (0, S(s)p) 

(o.3) VtX0 V(l < c(m, s) Av % p) * . 

Vxo KP) Vx KP) 

The following consequence of (0.3) is somewhat interesting. 

Let M m be a complete Kahler manifold with nonnegative holomorphic bisectional cur- 
vature. Suppose that M contains a compact subvariety V of complex dimension s. Then 
there exists C = C(m, s) > such that for S(s)p > p' ^> 1, 

\r I \ / \ 2(m— s) 



(0.4) < C 

VxoiP') \P 



In particular, 



lim -tt, r < oo. 

p^oo p 2 (m-s) 



The result sharpens the Bishop-Gromov volume comparison theorem in the presence of 
compact subvarieties. When M is simply-connected the result is in fact a consequence of 
the splitting theorem proved in [NT2, Theorem 0.4], via a completely different approach. 
It is not clear it can be derived out of any previously known result in the general case. For 
an entire analytic set V in C m (of dimension s) , one can define the Lelong number by 

n * r (V) (m - s) A v , X0 (p) 
^oo(V) = sup hmsup , z • 

x EM p^oo V x \P) 

Stoll showed that V is algebraic if and only if ^oo(V) < oo. The following result generalizes 
his result to analytic sets in curved manifold (we only succeeded in codimension one case). 

Let M be a complete Kahler manifold with nonnegative bisectional curvature. Let V be a 
analytic divisor of M . Then V is defined by a 'polynomial function' (holomorphic function 
of polynomial growth) if and only if VooiV) < oo. 
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In Section 4 we also obtained other results, one of which generalizes the classical tran- 
scendental Bezout estimate for codimension one analytic sets. (It has been known that the 
result fails for the high codimension case by the famous example of Cornalba and Shiff- 
man [CS], even for the Euclidean case.) We found the connection between the parabolic 
approach, especially the Li-Yau-Hamilton estimate, and the classical Nevanlinna theory 
interesting and believe in that the parabolic approach should be the most nature/effective 
approach in extending sharp results from Euclidean spaces (linear) to the curved complex 
manifolds (nonlinear, in certain sense). 

Finally, we discussed the relation between the Li-Yau-Hamilton inequality proved in this 
paper and the previous computations in [FIN] on the reduced distance and the reduced 
volume modelled on Ricci expanders. In particular, we prove a sharp lower bound on 
the heat kernel for the time dependent case, which, in a sense, is dual to Perelman's 
monotonicity of the reduced volume. We also explain how one can view this result, and 
more importantly, Perelman's monotonicity of the reduced volume as a nonlinear version 
of earlier work of Cheeger-Yau [CY] and Li-Yau [LY] on the heat kernel estimates for heat 
equations/Schrodinger equations. We also proved several local monotonicity formulae for 
Ricci flow on the forward reduced volume defined in [FIN] and Perelman's entropy, without 
any curvature sign assumption. This again follows the observation of Ecker in [Ec2] , where 
he derived a localized version of Huisken's monotonicity formula for the mean curvature 
flow. As an application, we prove that the blow-down limit of a so-called type III solution 
to Kahler-Ricci flow with bounded nonnegative bisectional curvature must be a gradient 
expanding soliton. This is, in certain sense, dual to Perelman's result on ancient solutions. 

Here is how we organize the paper. In Section 1 we prove the interpolating version of 
the new Li-Yau-Hamilton estimate, which in particular, includes (0.1). In Section 2 we 
derive some monotonicity formulae and the heat kernel comparison theorem. In Section 
3 we show the interpolation between the new Li-Yau-Hamilton inequality and Perelman's 
entropy formula in [P] . In Section 4 we derive the localized version and prove the manifold 
version of Stoll's theorem. In Section 5, we discuss the relation of the new inequality 
with the work of [FIN] , formulate a new matrix Li-Yau-Hamilton expression on curvature 
tensors for Ricci flow, and show a sharp heat kernel lower bound estimate. 

Acknowledgement. The special thanks go to Professor B. Chow as one can see that his 
contribution is crucial to a couple of results in this paper. He generously encouraged the 
author to publish the results alone even though it should really be a joint paper. The author 
would also like to thank Professors H.-D. Cao, Tom Ilmanen, Peter Li and Jiaping Wang 
for helpful discussions, Professors Klaus Ecker, Luen-Fai Tarn, H. Wu for their interests. 

§1 A new matrix Li-Yau-Hamilton inequality for Kahler-Ricci flow. 

Let M m be a complete Kahler manifold of complex dimension m. Let (M,g(t)) be a 
solution to Kahler-Ricci flow: 

d 

Here R a p{x,t) is the Ricci tensor of the metric g a p(x,t). Let u be a positive solution to 
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the forward conjugate heat equation: 



(1.2) 



( 



d_ 

di 



A J u(x, t) = 1Z(x, t)u(x, t). 



Here TZ(x, t) is the scalar curvature. We shall prove the following new matrix Li-Yau- 
Hamilton/differential Harnack inequality. 

Theorem 1.1. Let (M,g(t)) be a solution to (1.1) with nonnegative bisectional curvature. 
In the case that M is complete noncompact, assume further that the bisectional curvature 
is bounded. Let u be a positive solution to (1.2). Then 



for any (1,0) vector field V . 

The assumption on (M,g(t)) having bounded nonnegative bisectional curvature can be 
replaced by (M, g(0)) having nonnegative bisectional curvature and (M, g(t)) with bounded 
curvature, thanks to the result of Bando [B] , Mok [M2] and Shi [Sh2] . The same applies 
to other results of the paper. We state the result under the stronger assumption just for 
the simplicity. 

Recall that in [CN] the authors proved that for a fixed Kahler metric (M, g a p{x)) with 
the nonnegative bisectional curvature and the positive solution u to the heat equation 
( — A) u = one has the matrix Li-Yau-Hamilton inequality: 



Therefore, one can think (1.3) is the nonlinear version of (1.4). In fact we shall show that 
there exists a linear interpolation between these two inequalities. The similar interpolation 
was established by Ben Chow [Chi] originally for the Li-Yau's gradient estimates for the 
heat equation and Hamilton's differential Harnack for the Ricci flow in the case of surfaces. 
We shall show such interpolation between the matrix differential inequalities (1.3) and (1.4) 
for any dimensions. In [N4] we have shown another family of Li-Yau-Hamilton inequalities 
which also serves as a higher dimensional generalization of Chow's result. 

Let us first set up the notations. For any e > 0, we consider the Kahler-Ricci flow: 



(1.3) 



U ap + Tdap + V-R a p + U a Vp + UpV a + uV a Vp > 



(1.4) 



U ap + T9 a p + U a Vp + UpV a + uV a Vp > 



(1.5) 



d_ 
di 



Consider the positive solution u to the parabolic equation: 



(1.6) 




We shall call (1.6) forward conjugate heat equation. We shall prove the following general 
interpolation theorem. 
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Theorem 1.2. (Chow-Ni) Assume that (M,g(t)) has nonnegative bisectional curvature 
satisfying (1.5). In the case that M is complete noncompact, additionally assume that the 
bisectional curvature is bounded. Let u be a positive solution to (1.6). Then 

(1.7) + -# Q ^ + euR a0 + u a Vp + upV a + uV a Vp > 
for any (1, 0) vector field V . 

Remark 1.1. Since Ben Chow helped with the proof at a critical step it is well-justified to 
attribute the above theorem as a joint result. 

It is easy to see that Theorem 1.2, serving an interpolation between (1.3) and (1.4), 
implies Theorem 1.1 and the earlier result (1.4) for the linear heat equation. Therefore, we 
only need to prove Theorem 1.2. The proof consists of the following several lemmas. 

Lemma 1.1. 

(1.8) (Jtyap = A-R Q/ 3 + R a p^fiRfi^ - R a pR p p. 

Here -R Q/ 3 7< 5 is the bisectional curvature. 
Proof. The second Bianchi identity yields 

(1-9) l^afi — \R"lV,aj3 = Ray,yP 

Ray,py ~ RpyypRap ~\~ R a pypRpy 
= -^of/3,77 ~~ RppRap + R-ap-ypRpj- 

The commutator formula gives that 

(1-10) R a p^ = R a p^ — R p p^R a p + R a p^R p p 

= Rap,-yy ~ RppRcep + RapRpp- 

The lemma now follows from the definition Ai? a/ g = | (i? Qj g 77 + R a p ^ . 
Lemma 1.2. Let u(x,t) be a solution to (1.6). Then 

(1.11) (j^ - u a p = R a p^UyS - \ (R a pu p p + Rppu ap ) + e(TLu) a p. 

Proof. Differentiate (1.6) we have 

(1-12) K) 7 5 = Rp a y5 u aP + 9 ^^ocPyl + e(Ku)^. 

By definition Au a p = \ (u a p,yy + u ap,yy)i m normal coordinates at a point. We need to 
calculate the difference between the partial derivative u a p 1 g and the covariant derivative 
u a p 1 g. Direct computations show that, for normal coordinates at a point 

(1-13) u yS,aP = u y5ap + U ssR a pys- 
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Using the fact that 

(1.14) u j5,cece = U ~f5,aa + ^"yp U p8 ~ Rp5 U ~fP 



and (1.13) we have 



(1.15) A*U 7 5 = -(-U 7 5 )Q5 + « 7 5,aa) 



1 

1 

2 



Combining the above with (1.12), we conclude that satisfies (1.11). 

The direct calculations give the following lemma. 
Lemma 1.3. Let (M,g(t)) be a solution to (1.5). Letu(x,t) be a positive solution to (1.6). 



(1.16) 

V at 



\ /U a Ug\ (Ku) a Up (JZu)gU a ,U a Up U a sUg s +U as Ug- s 

A = e + e enZu) — ^ — 

J \ u J u u ir u 



-2 
+ 



it 2 -u 

U as UgUp + Up s U a Us + U a sUpU s + Up s U a U s 

2 ' 



it" 



(1-17) I — - A I (^-^ J (x, £) = e— ^ - - e-R a 0. 

(1.18) 

^— - A^j (euR a p) = e 2 HuR a p + eu (R a ^R^ s - R^pRpp) 

e (V s uV s R a p + V- s uV s R a p) 
+ (e 2 - e)w (R ce p J gR^s - R a pR p p + A-R a ^) . 

Proof. In the derivation of (1.16), a commutator formula has been used. In the derivation 
of (1.18) we have used 

d 

We also need the following result which is an easy consequence of Cao's differential 
Harnack for the Kahler-Ricci flow [CI]. 
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Lemma 1.4. Let (M, g(t)) be a complete solution to (1.5) with nonnegative bisectional cur- 
vature. In the case that M is complete noncompact we further assume that the bisectional 
curvature is bounded. Letu(x,t) be a positive solution to (1.6). Then 

(1.19) AR a p + R a fafiRy S - (—^—VsR a p + -^—V 8 R a p \ + R^i— ^— — 5 — + > 0. 

V C Uj C Uj J C U C Hi c L 

Once we have Lemma 1.1-1.4 we can give the proof of Theorem 1.2. 
Proof of Theorem 1.2. Let 

^ V a/3 — U a/3 + j9a/3 ~ 

and 

N a p = N a p + euR a p. 

By taking the minimizing vector in (1.7), one can see that Theorem 1.2 is equivalent to 
N a /3 — 0- Since the maximum principle needs some growth conditions for the noncompact 
manifolds, we first prove the theorem for the case that M is compact. 

Compact case. By Lemma 1.1-1.3, we have that 
(1.20) 

~ = R ap~<y5 U l5 ~ ^ { R otpU p p + R p pU ap ) + t{Ku) a p 

- e e h e(Ku) — ~^ H — 

u u u z u 

^ u a u p\ u s\ 2 1 RasU s Up + R s pUsU a 

U 2 U 

U as 1lgUp + Up s ll a Us + U a sUpU s + U0gU a U s 

u 2 

IZu u u 

+ e T^ uR a p + eu ( y R a p^Rjs — R apR p p) 

- e (V s uV- s R a p + V- s uV s R a p) 

= R ap 1 s^i5 - 2 (RapN p p + R p pN af ^ + -N ap N p p - jN a p 

If U a U p \ ( UpUp\ £ 2d d 

+ - ( u ap J I u p p J + alN^p - ue R a pR p p 

Ui \ Uj / \ Uj ' 



+ e 2 uY r 



a(3 



where 



V - AD 7? D f^ sU Y7 Z? ^ 5W V7 Z? \ D VjUVsU R a 
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By Lemma 1.4 and (1.20) we have that 
(1.21) 



§- f ~ a) N a (3 = RafhsNyS - \ (R a pN p - p + R p pN ap ) + (ell - ^ iV a/3 



+ ^ N ap { N pf3 ~ £ U Rpf3) + ^ (Nap - e-ul^ap) iV^g 
1 / U a U p \ / 



> R a fagN^ s - - (R ap N p p + R p pN ap ^j + ^eTZ - iV a ^ 

+ 77- ^ap - CU-R^g) + 77- (Nap - e^ap) N p p. 



Using the observation that the right hand side of (1.21) satisfies the null-vector condition of 
the tensor maximum principle of Hamilton [H2], we have proved the case M being compact. 

Noncompact case. First recall the fundamental derivative estimate of Shi. For g a p(x,t), 
a solution to (1.5) on M x [0, T] with bounded (in space-time) nonnegative bisectional 
curvature, there exist A k > such that 

(1-22) l|V fc ^ 75 -|| 2 < ^ 

on M x [0, T] . For our consideration we only need to prove it for T small. The estimate 

(1.22) is proved in [Shi]. For the sake of simplicity, we will show the matrix differential 
Harnack inequality Theorem 1.2 for the case e = 1 under the above assumption (1.22). In 
fact, what needed is (1.22) for k < 2. We also need the perturbation trick from [NT1] (see 
also [CN]). Namely we first shift t by 26, where 6 is a small positive number. After the 
shifting we can have estimates on *u, |Vit| and \u a p\. The goal is to show that there exits 
a b > such that 

(1.23) / / exp(-6(r^(x) + 1)) (- + + K«| 2 ^) dfidt < oo. 
Js Jm \u u J 

Here tq(x) is the distance to x from a fixed point o G M with respect to the initial metric. 
We need the estimate (1.23) to apply the tensor maximum principle from [N5] (see also 
[NT2] for the original time-independent version). 

In order to get control on u (or ^) first we need the following Harnack inequality of 
Guenther [Gu]. 

Theorem 1.3. (Guenther) Let (M,g(t)) be a solution to Ricci flow satisfying (1.22). Let 
u be a positive solution to the forward conugate heat equation ( — A) u = TZu. Then for 
sufficient small T (only depending on A k ) there exist a, B\ > only depending on A k such 
that 

(1.24) u{ Xl M) < u{x 2 ,t 2 ) (^) 2 exp( r -^^^- + B 1 (t 2 -t 1 )). 

\tij ct(t 2 -ti) 
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for any T > t 2 > t\ > 0. Here r(xi,x 2 ,ti) denotes the the distance between x\ and x 2 
with respect to the metric at t = t\. 

The result above was proved in [Gu] through a gradient estimate of Li-Yau type on 
compact manifold. Since one can apply the localization techniques as in [LY, page 161] 
(see also [NT1, page 647], [Shi] ) one can easily generalize the gradient estimate, thus the 
above Harnack estimate, proved in [Gu] to complete noncompact case. From (1.24), one 
can deduce that for small 5 there exists a constant b 2 ,B 2 > 0, where b 2 = b 2 (A]~, 5) and 
B 2 = B 2 {u(o,%),u(o,T- %),A k ,6) such that 

(1.25) (J- + i?J {x, t) < B 2 exp(6 2 (r2(x) + 1)) 

for (x,t) G M x [|,T- §]. 

Observe that we have the following two equations. 

(1.26) (!__A) u 2 = Ku 2 -\Vu\ 2 



dt 



and 



(1.27) (J^ - A^ \Vu\ 2 = -K^| 2 - \u a0 \ 2 + (V(%«),V«) + (Vu, V{Ku)). 

The desired estimate (1.23) follows from (1.26), (1.27) by the argument through integration 
by parts in [CN, Lemma 3.1]. Once we have established the estimate (1.23) one can apply 
the perturbation argument as in [NT1] (see also [CN]) together with the tensor maximum 
principle in [NT2, Theorem 2.1] or [N5, Theorem 2.1] to conclude the proof of the matrix 
Li-Yau-Hamilton estimate (1.7) for the complete noncompact case. Note that Theorem 
1.3 only provide the estimate for short time. But we can iterate the argument to prove 
the result for all time. An alternative is that once one has the upper bound estimate 
at some earlier time one can also make use of the heat kernel estimate in [N5] for the 
time dependent heat operator (and uniqueness of the positive solution) to get estimates of 
the positive solution for the later time. We also should remark that the matrix Li-Yau- 
Hamilton (1.7) gives a sharp Harnack (which is more precise, compared with Theorem 1.3). 
Please see Corollary 2.1 in the next section. However, we do need the rough estimate in 
the proof to apply the tensor maximum principle. 

Corollary 1.1. Let u{x,t) be a positive solution to (1.6). Then 

ill 

(1.28) Alogw + eft+ — > 0. 

t 

If the equality holds for some (xo,to) with to > 0, then (M,g(t)) is an expanding Kahler- 
Ricci soliton for the case e > 0, and (M, g) is isometric to C m for the case e = 0. 

Proof. Let 

ill 

Q = Alogw + eft+ — . 
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and Q = g a ^Y a p. Since N a p > and Q = g a ^N a p we have that N^(x ,t ) = 0. By the 

strong maximum principle we know that N a p = for all t < to- N a p = is nothing but 
the equation in the definition of the gradient expanding soliton. 

For the case e = 0, we apply the same line of argument. In this case we have Q = 
since N a p = N a p = 0, for t < to- Now from the equation 

= (| - A) eQ = tiR^loguUlogu)? + t (u ap - ^) [u m - 

It 

> o 

we have that (\ogu) a p = too. From the definition of iV a/ 3 we then have 
and 

(logw) Q/3 = 

from which it is easy to see that (M, g) is flat and in fact isometric to C m , since that the 
curvature (see, for example, [KM, page 117]) can be written as 

R - d 4 (*/) , a «,( d*(tf) \ f d\tf) 

a ^ s dz a dz~Pdz^dz~ 5 \dz*dz a dz~< ) \dzPdz~Pdz~ 5 

and that / = log it (defined to be) is a convex function. 

Remark 1.2. The case e = case in the above Corollary 1.1 is just the original Li-Yau's 
estimate [LY], which holds with nonnegativity of the Ricci curvature. The equality case 
implies the manifold is lR n is implicit in the proofs of [LY] and proved explicitly in [N3] for 
Riemannian manifolds with nonnegative Ricci curvature. It would be interesting to see if 
Corollary 1.1 is true, for e > case, for complete solutions to the Kdhler-Ricci/ Ricci flow 
without assumptions on the sign of the curvature. 

When the manifold is compact, since we known that, by passing to its universal cover, it 
is products of C k with compact Hermitian symmetric spaces. Without the loss of generality 
we can assume that the first Chern class c±(M) is a positive multiple of the Kahler class. 
Then the Kahler- Ricci flow will have singularity, say at t = 1, and the normalized flow has 
long time existence. The rescaling is given by g = yz^-g an d the reparametrization is given 
by s = — log(l — t). Therefore, Theorem 1.2 has the following equivalent form. 

Theorem 1.2'. Let M be a compact Kahler manifold as above. Let g(x,t) be a solution 
to (1.1), and let u(x,t) be a positive solution to (1.2). Assume that g(x,s) is the solution 
to the normalized Kahler Ricci flow. Then 

(1-29) (logu) a ,g + R a p + ^—[9af3 > 0- 

Here R a p is the Ricci tensor of g. 
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§2 Monotonicity formulae. 

In this section we derive some monotonicity formulae out of the matrix Li-Yau-Hamilton 
estimate, as well as its trace, proved in last section. In order to make the argument unified 
for both cases, with and without Ricci flow, we work with the interpolation version (1.7). 
Let (M,g(t)) be the solution to the Ricci flow (1.5) and let u(x,t) be the positive solution 
to (1.6). 

The first result is the monotonicity of the partition function (also called Nash's entropy 
in [FIN]) defined by 

(2.1) J\f(g,u,t) = — / ulogudv — mlog(nt) — m. 

J M 

Simple computation shows that 

(2.2) ^T~= f (-Alogu-eft- —)udfH < 0. 

Here dji t is the volume element of g(t). The following result is a direct consequence of 
Corollary 1.1. 

Proposition 2.1. Let (M,g(t)) be a solution to (1.5) with nonnegative bisectional curva- 
ture. Let H(x,y,t) be the fundamental solution to (1.6). Then 

±Af(g,U,t)<0 

for any positive solution u to (1.6) and 

— / H log H dfi t — mlog(nt) < m 
Jm 

with the equality holds for some positive t if and only if the manifold is an expanding 
gradient soliton (isometric to C m in the case e = 0). 

For the fixed Riemannian metric case the above was proved earlier in [N2] for the 
manifold with only nonnegative Ricci curvature. We believe that the result should hold for 
Ricci flow even without assumptions. But at this moment we can only prove it for Kahler- 
Ricci flow through the matrix Li-Yau-Hamilton inequality, which assumes the nonnegativity 
of the bisectional curvature. The result above gives a characterization of expanding solitons 
using the partition function. The similar formulation for the shrinking solitons also works 
for the partition functions related to Perelman's entropy formula [P]. Namely, consider the 
backward Ricci flow -§^gij = 2Rij on M x [0, To], where M is a Riemannian manifold of 
real dimension n. Let u(x, r) be a solution to the backward adjoint heat equation 

(1.6') (-^-A + n\u(x,r) = 0. 



14 



LEI NI 



Similarly one can define 

~ f Ti TI 

(2.1') Af(g,u,r) = - ulogud^r - - log(47rr) - -. 

Jm 1 1 

In Proposition 1.2 of [P], Perelman proved that 

(2.2-) f =/ M (-A.o gM+K -£)^<0. 

The dual version of Proposition 2.1 states as follows. 

Proposition 2.1'. Let r) 6e £/ie fundamental solution to the adjoint heat equation. 

Then 

f TI TI 

- / H log Hdfi T - -log(47rr) < - 

with equality holds (or in (2.2') ) for some positive r if and only if (M, <7y(r)) is a gradient 
shrinking soliton. 

Proof. Since we do not have Corollary 1.1 in this case we need other arguments. In fact, 
tracing the equality case of the proof of Proposition 1.2 of [P] we have that Rij—ViVj log H 
is diagonal. Namely Rij — V^Vj \ogH = n ~ A ^°s H g-. On the other hand, the equality in 
(2.2') further implies that — V^Vj logH = -^gij. 

In [N2], the relation between the value of J\f(t) (as well as the entropy functional W) 
as t — > oo and the asymptotic volume ratio at infinity (also called the 'cone angle') was 
proved for the linear heat equation. The similar relation, between the ^-constant in the 
K-non-collapsing of volume defined in [P] and the large time limit of the partition function 
Af(g, u, r) (as well he entropy functional W) should also be true for ancient solutions to 
Ricci flow. 

Another application of Theorem 1.2 is an entropy monotonicity for the ancient solutions. 
Let (M, g{t)) be an ancient solution to (1.5) and let uhe a positive solution to (1.6), both 
defined on M x (— oo,0]. Theorem 1.2 implies that 

( l ogu) a p + R a p > 0. 

From this one can easily see that 



Af(g,u,t) := -/ ulogu dfi t 
Jm 



is monotone non-increasing. 

If M is compact, one can obtain such u by taking limit of solutions with initial data at 
t = — i as in [FIN] (where the immortal solution is studied). More precisely, let Ui(x,t) be 
a solution to (1.6) with u(x, —i) = yjpij ■ Letting % — > oo, one can extract a limit u^ > 
(since f M Ui = 1 and Ui > 0) which is defined on (— oo, 0]. In this case fj^Uoodv = 1. 
Applying the Jensen's inequality one has that 

A%, w^, t) < logV (t). 



LI-YAU-HAMILTON INEQUALITY 



15 



The right hand side is another monotone non-increasing quantity along the flow. When M 
is complete noncompact one can obtain such a u similarly by solving Ui(x,t) with initial 
condition at t = — i and anchoring Ui(o, — 1) = 1, where o G M is a fixed point. In fact one 
can even get Ui integrable by taking it to be scalar multiple of the fundamental solution 
(with initial data being the delta function at t = —i). 

The second application is to derive the heat kernel comparison theorem and Huisken 
type [Hu, El-2] monotonicity formula for the analytic subvarieties in M. We start with 
the case e = since the results seem to be more useful at this moment. 

Theorem 2.1. Let M be a complete Kdhler manifold with nonnegative bisectional cur- 
vature. Let H(x,y,t) be the fundamental solution of the heat equation. Let V C M be 
any complex subvariety of dimension s. Let Ky(x, y, t) be the fundamental solution of heat 
equation on V . Then 

(i) 

(2.3) K v (x,y,t) < {7rt) m - s H{x,y,t), for any x,yeV. 

If the equality holds, then V is totally geodesic. Furthermore if M is the universal 
cover of M with covering map n and V = 7r _1 (V), then M = Mi x C k for some 
Kdhler manifold Mi which does not contain any Euclidean factors, with k > m — s. 
Moreover V = Mi x C l with I <k. 

(ii) 

(2.4) - J^t) m - s H{x,y,t)dA v {y) > 0, for any x G M. 

Similarly, if the equality holds for some x G M at some positive time t, then 
M = Mi x C k with k>m-s. 

Proof. For any smooth point y G V, choose a complex coordinate {z\ ■ ■ ■ , z m ) such that 
(z\,--- ,z 8 ) is the coordinates for V. Let i, j, k,--- denote the coordinate on V and 
a, 6, c, • • • denote the coordinates in the normal directions. Then we compute 

(2.5) 

(A^ - {(ntr~ s H(x, y, t)) = (Am - g al V a V~ h - ^ {t m ~ s H(x, y, t)) 

= (nt) m - s U M H -H t - g al V a V- h H - 



-{TTt) l-VaV^iiH — -Hg ah \g 

- W 



iV^tfl 2 



By Theorem 1.2 we have that 



4 W) -^) {(*t) m -'H(x,y,t))<0. 
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Noticing that for x G V, \Sm. t ^a( , Kt) rn ~ 8 H(x, y, t) |y = S x (y). This proves that (2.3) by the 
maximum principle, (2.4) by the integrating (2.5) on V. (For the case V is singular, one 
can refer to [LT] for the justification on the validity of the integration by parts.) 

In the case when the equality holds in (2.3) we have that (7Tt) m ~ s H(x, y, t) satisfies the 
heat equation. Hence equality holds in (2.5) for any x, y G V. This implies that 

(2.6) |V^ log#| 2 = 0, for x, y G V, t > 0, 

which implies that < Vr 2 (i, y), v >= for any x, y G V, and any normal direction v. Here 
we have used the fact that lim t ^ —t log H(x, y, t) = r 2 {x,y). See [CLY1], for example. 
This implies that for any smooth point x G V, any minimizing geodesic starting from x lies 
totally inside V. More precisely, let r ) v {s) be a short minimizing geodesic emitting from x 
with 7^,(0) = v such that v G T X V. Denote by h be the distance function from V. Then 

±-h{ lv {s)) = 2Re(VK 1 ' v {s)) 

= -2Re(Vh, S y v (s)) 
s 

= -Re{Vh,Vr 2 {x^ v {s))) = 0. 
s 

Here V/i = V a /ig|^ and j' v (s) = dz d ^ gf^- Therefore, V is totally geodesic. This in 
particular shows that there is no singular points in V, which rules out the possibility that 
V is union of several totally geodesic submanifolds with singular intersections (in which 
case the heat kernel comparison (2.3) has strictly inequality). Moreover, by lifting the 
computation to the universal cover we can assume that M is simply-connected. Then the 
manifold M splits by Theorem 0.1 of [NT2], more precisely, Theorem 2.1 and Corollary 

2.1 of [NT2], since N a p = VctV pH H h jHg a p > and its null space is at least 

of m — s dimension, by (2.5), for any x, y G V and t > 0. Notice that N a p does not 
satisfies the linear Lichnerowicz heat equation. However the inequality (1.21), satisfied by 
N a p (since e = 0), is enough for the argument in the proof of Corollary 2.1 of [NT2]. Then 
M = Mi x M2 such that the tangent space of M2 consists of the null space of N a p. The 
factor Mi is isometric to C m_s follows from the same argument in the proof of Corollary 
1.1, since — V a Vg-f/"+ VaH ^b H _ ±Hg a i = 0. (One can also use Corollary 1.3 of [N2, page 
331].) More precisely, if we write a point x G M as x = (x±, x 2 ) according to the splitting, 
we can write the heat kernel H(x,y,t) = Hi(xi,y\,t)H2(x2,y2,t). Then on M 2 we have 
that (logH^O^) J/2)) a /3 + jdap = by the definition of the splitting. Therefore one can 
apply Corollary 1.1 to conclude that M 2 = C fc . If (2.4) holds equality for some x G M, it 
implies that the right hand side of (2.5) is zero. Then the argument as above also applies. 
Note that we may not have V totally geodesic since x may be a singular point. 

Remark 2.1. 1) In the case V is not smooth, the existence on K(x,y,t) was justified in 
the work of Li and Tian [LT]. They also obtained a similar upper bound estimate as (2.3) 
for the special case M = P m with the Fubini-Study metric using very different method. 
Their method produces better upper bound for the special case M = P m . Our estimate here 
works for general Kahler manifolds with nonnegative bisectional curvature. 
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2) The similar heat kernel comparison was first proved in [CLY2] for minimal subman- 
ifolds in space forms by quite different method. The result in Theorem 2. 1 is more general 
than [CLY2] in the sense that it holds for any Kahler manifolds with nonnegative sectional 
curvature in stead of space forms. However it is also more restrictive since it only applies 
to analytic subvarieties. 

3) In the part (ii) of Theorem 2.1, the manifold M\ may not be V. This could happen, 
for example, in the case M = C m and V is a union of two hyper-planes and x lies on 
the intersection subvariety. If we further assume that equality holds for all smooth points 
x G V, we do have the same conclusion as part (i). 

4) The monotonicity (2.4) is enough for applications in [N4]. Namely, one can prove the 
comparison results on the dimensions of polynomial growth holomorphic function spaces ob- 
tained in [N4], using (2.4) in stead of the other Li- Yau- Hamilton inequality proved therein. 
In fact one can derive the results in [N4] through the following corollary, which is a special 
case of Theorem 2.1 (or Corollary 2.2) of [N4]. It is however enough for the applications 
considered in [N4] . 

Corollary 2.1. Let f be a holomorphic function. Let V = Z(f), the zero locus of f. 
Denote 

(2.7) w(x,t) = [ H(x,y,t)Alog\f\ 2 (y)dfi(y). 

JM 

Then 

(2.8) tw(x, t) = {Tit) [ H(x, y, t) dA v (y). 

Jv 

Moreover 

(2.9) ^(tw(x,t))>0. 

If the equality holds for some point x G M and some positive time t, then the universal 
cover (of M) M splits at least a factor of C 

Proof. Notice that A = g af3 dz f dzl3 , which differs a factor of 4 from [N4]. Let / be a 
holomorphic function defined on M. Here we do require that / does not growth too fast. 
For example, it will be sufficient if / is of polynomial growth or is of finite order in the 
sense of Hadamard (see (3.1) of [N4] for precise definition). We can write 

v(x,t) = [ H(x,y,t)log\f\ 2 {y)d f i(y) 

J M 

as a solution to the heat equation with initial value log \f\ 2 (y). The requirement on / is 
to make such representation formula of v(x,t) meaningful. Then w(x,t) is defined to be 
■§iv(x, t) as in [N4]. By the definition of w(x, t) in (2.7) it is easy to see that w(x, t) is also 
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a solution to the heat equation with initial data given by the measure A log \f\ 2 . It is easy 
to see that w(x, t) = -§j;v(x, t). By the Poincare-Lelong formula we know that 

(Trt) H(x, y, t) {^dd log \f\\y)^ A ^_ j ); = (nt) H(x, y, t) dA v (y). 

On the other hand, the direct calculation shows that the left hand side of the above equation 
is equal to 

tf H(x,y,t)Alog\f\ 2 (y)^=t [ H(x,y,t)Alog\f\ 2 (y) dfi(y) 

JM m - JM 

= tw(x, t). 

This proves the first statement of the corollary. The monotonicity (2.9) is just a special 
case (ii) of Theorem 2.1. The proof above also shows that w(x,t) has the same meaning 
as in Lemma 3.1 and Theorem 3.1, Theorem 4.1 of [N4]. 

The case with Ricci flow, namely e = 1, can be formulated similarly. In order to do so 
we have to explain some notations. For a time-dependent metrics deformed by the Kahler- 
Ricci flow equation (1.1), we call H{x, y, t, t ) a fundamental solution to the heat equation, 
if for any u(x, to) = <p(x), where <p(x) is a compact-supported smooth function, the solution 
u(x,t) to the heat equation (^ — A) u(x,t) = with initial condition u(x,to) = (p(x) is 
given by the formula 

u(x,t)= H(x,y,t,t )(p(y)dfj,(y). 

JM 

It is easy to check that H{x, y, t, to) must satisfies the forward conjugate heat equation (1.2). 
When the meaning is clear in the context (mostly to = 0) we just simply write H(x, y, t, to) 
as H(x,y,t). Adapting the notation used in the proof of Theorem 2.1, the restriction of 
the Kahler-Ricci flow (1.5) on V reads as ^gij(x,t) = —Rfj(x,t). We denoted g lJ ' by 
TZ\>. It is easy to see that 1Z\> is well-defined. Therefore the fundamental solution to the 
heat equation on V, K v (x 7 y,t) satisfies 

(2.10) ^ - v(x, t) = K v {x, t)v{x, t). 

We have the following comparison and monotonicity result. 

Theorem 2.1'. Let M be a complete Kahler manifold with bounded nonnegative bisectional 
curvature. Let H(x,y,t) be a fundamental solution to the heat equation on M. Let V be 
a complex subvariety of M of dimension s. Let Ky(x,y,t) be the fundamental solution to 
the heat equation (with respect to the induced metrics) on V. Then we have (2.3) and (2.4). 
Moreover, the equality (for positive t), in either cases, implies that the universal cover (of 
M) M has the splitting M = M x x E k , where E k is an gradient expanding Kahler-Ricci 
soliton of dimension k > m — s. 
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Remark 2.2. One can think (2.4) as a dual version of Perelman's monotonicity of the re- 
duced volume since the reduced volume of [P, Section 7] is, in a sense, a'weighted volume' 
of M (with weight being the heat kernel of a 'potentially infinity dimensional manifold' 
restricted to M , as explained in Section 6 of [P]), while here the monotonicity is on the 
'weighted volume' of complex submanifolds with weight being the heat kernel of M restricted 
to the submanifold. The reduced volume monotonicity of Perelman has important appli- 
cations in the study of Ricci flow. We expect that (2.4) will have some applications in 
understanding the Kahler-Ricci flow and its effect on the complex geometry of analytic 
subvarieties. 

Taking the trace of the matrix estimate in Theorem 1.2 and integrating along the space- 
time path as in [LY] , we can have the following Harnack estimates for the positive solutions 
to the forward conjugate heat equation. This gives a sharp version of the previous rough 
estimate of Guenther in [Gu]. 

Corollary 2.1. Let (M,g(t)) be a solution to Ricci flow (1.1) and u(x,i) be a positive 
solution to (1.2). Then 

|V^| 2 _ ut rn > Q 
u 2 u t ~ 

and for any ti > t\, 

(2.11) u{x 2 ,t 2 )t 2 n > u(xi,*i)ty*exp (~ mf ^ 2 

Here ^(t) is a path with 7(^1) = x\ and 7(^2) = x 2 . 

Note that we do not have the factor 4 due to our choice of A and the gradient |V/| 2 = 
g^f a fp. Also 

. / , , .0 dz dz^ 

IVWI^M-s-s- 

We list this consequence here since it implies the monotonicity of t m u(x,t). 

§3 Interpolation between Perelman's entropy formula and the new Li-Yau- 
Hamilton inequality. 

The purpose of this section is two-folded. First we give a different proof of Theorem 
1.2. The second purpose is to show that the by-product of this different proof also implies 
Perelman's monotonicity of entropy, as well as the energy. The computation in this section 
has its real version. See [Ch2] and the forth-coming book [CLN]. The main computation 
is summarized in equation (3.7) below, which we call a pre-Li-Yau-Hamilton equality. The 
equation (3.7) can also be viewed as a matrix version of Perelman's entropy monotonicity 
formula. In a sense, one can think that the matrix Li-Yau-Hamilton inequality proved in 
Section 1 is dual to Perelman's entropy formula in [P, Section 3]. 

Consider the Kahler-Ricci flow: 



8 

(3-1) = e "^W 
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where e is a parameter and the conjugate heat equation: 
(3.2) (J^-- A + eUj u(x,t) = 0. 

When e < 0, (3.1) is a forward Ricci flow equation and (3.2) become forward conjugate 
heat equation. The equations look different from those in Section 1 since in this section the 
case of e < corresponds to the forward Ricci flow and the case of e > corresponds to 
the backward Ricci flow. For example e = 1 is exactly the setting for Perelman's entropy 
and energy monotonicity. Notice that (3.2) becomes the backward conjugate heat equation 
for e = 1. For the positive solution u(x, r) we define the (1, 1) tensor Z a p by 

Let Ax, denote the Lichnerowicz Laplacian on (1, 1) tensors, which is defined by 

for any Hermitian symmetric (1.1) tensor r} a p. It is known, see for example [CI], that 

The direct calculation as in [NT1, Lemma 2.1] shows that 
Lemma 3.1. For any C 2 -function f(x, r) 



(3-3) ( £ " ^ ) /«? 



a/3 



Remark 3.1. The similar result as the above lemma hold for Ricci flow on Riemannian 
manifolds. Please see [CLN] for details. 

Now with the help of Lemma 3.1 and 1.1 we can calculate (^-Ax) Z a ^ as follows 
using the equation ( — Ax,) (log it) = — eR, + |Vlog"' 12 



u 



(3.4) 

= e(K) a0 - (^(log U ) 7 (log«)j) -e(l + e)A L R a0 

= -R a p 7 5( lo Su)j(logu) s - (logw) Q7 (logw)^ - (logw) Q7 (logw) 7 ^ 
- [Qogu) a p] (logu) 7 - [(\ogu) a p\ (logu) 7 - t 2 A L R a - p . 
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Hence 
(3.5) 

(Jff ~ A l) Zap = -e 2 (&R a j5 + Rajihs&ys + V 7 # a ^V 7 log u) + V 7J R Q/5 (-3-V 7 log u) 

+ ^ 7 5(^V^logw)(^V 5 logw)^ 

+ e 2 R a ^R^p - (logu) a7 (logu)^ - (logw) a7 (logw) 7 ^ 
+ V 7 (Z a/5 ) V 7 log u + V 7 (Z a/5 ) V 7 log u. 

Regrouping terms yields 



(3.6) 

(' 



^ - A L ^J Z a - p = -e 2 (aR o0 + R aBl ~ s R^ + V 7 # a ^V 7 log u) + V 7 # Q/5 (^V 7 logu) 

+ ^a^(^V 7 l0gw)(^V 5 l0gw)^ 



- (logu) a7 (logu)^ + V 7 (Z Qj g)V 7 logw + V 7 (Z a/5 )V 7 logw 
+ ^Z a7 (eR 1 p + (logu) 7( g) + - (ei? a7 + (logu) a7 ) Z^. 



Let 



and 



Z a p — Za(3 ~ ~9aP 



Y ap = A Rap + RapjSRjS +^ 1 R a p(^V^ log u) + V^R a p ( - V 7 log u) 
+ RaPjS ( - V 7 log U) ( - V S log U) . 

Notice that Y a p, defined after (1.20) in Section 1, is related to Y a p above through the 

~ Ft 

equation Y a p = Y a p (remember that — e here corresponding to e in Section 1). From 

(3.6), we can derive the equation for Z a p as follows. 

Lemma 3.2. (Chow-Ni) 
(3.7) 

( J~ ~ Al ) ^ap- = (Jfp - A i^J Z a - p + ^g a p - ^eR^ 

= ~ (f Y aP - ^ R ap) ~ (l0gw)a 7 (l0gw)^ 

+ V 7 (Z a ^)V 7 logw + V 7 (Z a ^)V 7 logw 

+ ^Z Q7 (eR 1 p- + (logu) 7/5 ) + - (e# a7 + (logu) a7 ) Z^ - ^Z a p. 
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Notice that Z a p = ^N a ^. With some labor one can check that (1.21) and (3.7) are 
equivalent. Namely one can derive one from the other, keeping in mind that — e here 
corresponds e in Section 1. One can also write (3.7) as 

( 3 - 7 ') (j; - z«e =(§p- Al ) z ^ + ~ \ tR ^ 

+ V 7 (Z Q/ j) V 7 log u + V 7 (Z a$ ) V 7 log u 

+ 7, + (logw) a7 - ^ga^j Z^. 

For e < 0, applying Lemma 1.4, the fundamental result of H.-D. Cao, we know that 
e 2 Y a p — ^R a p > under the assumption that M is a complete Kahler manifold with 
bounded nonnegative holomorphic bisectional curvature. Hence the tensor maximum prin- 
ciple and (3.7) implies that Z a p < 0, which is equivalent to the statement of Theorem 1.2. 
Namely, (3.7) does lead to another proof of Theorem 2.1. 

Even though the computation (3.7) and (1.21) are essentially equivalent, (3.7) has the 
advantage that when e = 1 it also implies Perelman's energy /entropy monotonicity for- 
mulae. The following is a more detailed computation of this claim. Let / = — logw, 
Z = g a ^Z a p. The tracing of (3.6) gives 



(3.8) 



(J^-A)Z = -eRspZrf - e 2 g^Y a - p - (/) Q7 (/) 7 * - V 7 ZV 7 / - V 7 ZV 7 / 



and 



g^Y a p = ATZ + R a pRap - V 7 ^V 7 /) - (^V 7 /)V 7 ^ + R a ^-J & ){-Jp). 

The following observation of Ben Chow is also useful. 
Lemma 3.3. (Chow) In the case e = 1, we have that 

(3.9) / (g^Yrf) udp= f (R a - (R^ + Up)) u dp. 

Jm v 7 Jm 

Proof. Follows from the integration by parts and the second Bianchi identity 1Z 1 = -R 7 q, q . 
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Remark 3.2. Please refer to [Ch2] and [CLN] for the Riemannian version of the above 
identity. 

Recall the definition of energy T . 

f (\Vu\ 2 \ 

Tig, u,t) = / h IZu da. 

Jm\u J 

Then (3.8) (with e = 1) and Lemma 3.3 implies the the following result. 
Proposition 3.1. (Perelman) 

(3-10) ^T( 9 ,U,t) = - f (|i^ + / Q/5 | 2 + |/ a/3 | 2 )^. 

Note that this is nothing but the energy monotonicity formula of Perelman in Section 1 
of [P]. We show that it follows from (3.8). 

Proof of Proposition 3.1. Let e = 1 in (3.8) we have that 

(3-8') (Jt - A^jZ = -g^Y a - p - (/) Q7 (/)^ - V 7 ZV 7 / - V^ZV^f - Z a pf p& . 

Now the result follows from direct computation on £pT(g, u,r) = £p J M Zu dp,, by applying 
Lemma 3.3. 

Similarly, if we trace (3.7) and denote Z = g af3 Z a p, we have that 
(3.H) 

(j^-A^Z= -eR^Z a - p - e 2 g^Y^ +"-11- (/) Q7 (/ ) 7 « - V 7 ZV 7 / - V 7 ZV 7 / 

+ ZafiitRpu - //3a) Z. 

T 

For e = 1, integration by parts as before gives 

(3.12) ^- [ Zudfj, T = - [ (\Z a p\ 2 + \f a p\ 2 ) udu T - - [ Zudfj, T . 

dr J M Jm v j t Jm 

The above equation is equivalent to Perelman's entropy monotonicity formula due to the 
following consideration. Let 

Af(g,u,r) := — / u log u d\i T — m log(7rr) — m. 
Jm 

Then Perelman's entropy 

W(g,u,r)= [ [r(2Af-\Vf\ 2 + n + f-2m\udfi T , 
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where f = — logu — mlog(7rr), can be expressed as 

(3.13) W(g,u,r) = ^-{tN) = t I Zud/j T +Af. 

dr J M 

Therefore 

(3.14) — W = r— / Zudn T + 2 Zud[i T 

dr dr J M Jm 



~ r / (\ z a p\ 2 + IM 2 ) udfi T 
Jm v 7 



which is nothing but the entropy formula of Perelman in [P, Section 3] . 

As pointed out in [P], there exists a statistical mechanics analogy of Perelman's entropy. 
If we identify the quantities above with the notation of [Ev2, Chapter I and VII], r is the 
temperature; —Af defined above is the log of the distribution function in [Ev2]; — W is 
the entropy S in [Ev2]; rj\f is the free energy F; is the energy E in [Ev2] (which is 

nonnegative by Proposition 1.2 of [P]) and the first equation of (3.13) is just the well-known 
equation S = — from thermodynamics. The negation of the right hand side of (3.14), 
measuring the deviation from an shrinking soliton, is called the heat capacity in [Ev2] . The 
entropy formula (3.14) implies the concavity of S in E, one of the defining properties for 
the entropy. This analogy also holds for the solution to the heat equation with respect to 
a fixed metric Riemannian metric with nonnegative Ricci curvature [N3] . 

Remark 3.3. For the case of e = 0, the similar computation as above gives the mono- 
tonicity formula in [N3]. The strange thing is that one can not get nice monotonicity 
formula in the case e ^ -1 or 0. Namely, the interpolation formula (3.6)/(3.7) does not 
give nice interpolation for energy / entropy after integration on M . One can view (3.6)/(3.7) 
as a matrix version of the energy /entropy monotonicity formula for e > 0. This partially 
answers one of the questions raised in the end of [N3] (still not satisfactory though). On 
the other hand, when e < 0, one can not get entropy monotonicity out of (3.7). Instead 
we have a pointwise Li- Yau- Hamilton inequality. However, for e = 0, both entropy and the 
differential Harnack follows from (3.6) (See [N3]J. The above discussion indicates that our 
new matrix Li-Yau-Hamilton inequality is dual to the entropy monotonicity of Perelman 
in some sense and there may perhaps be certain profound duality behind the scene. 

Another puzzling point is that so far we have not been able to verify a matrix Li-Yau- 
Hamilton inequality analogue to Theorem 1.2 for the Ricci flow on Riemannian manifolds, 
even though the above computation (3.7) holds for e > for the Ricci flow on Riemannian 
manifolds (which was done in [Ch.2]). In short, the interpolation between positive and 
negative e by now only works in Kahler category. The validity of the Riemannian case is 
pending on the verification of a new matrix Li- Yau-Hamilton estimate similar to Hamilton's 
famous work [HI]. Please see Remark 5.2 for further details. 

§4 A local monotonicity formula and its applications. 

This section is inspired by the work of Ecker in [El-2]. Let M be a complete Kahler 
manifold with nonnegative bisectional curvature (unless specified otherwise). In this section 
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we study the localization of the previous established monotonicity (in Section 2) for a fixed 
Kahler metric. (We leave the Kahler-Ricci flow case to a later discussion.) In [El], the 
localized monotonicity formula is proved for mean curvature flow in Euclidean spaces. Since 
we are dealing with curved spaces here, we need some extra ingredients, which includes 
Theorem 2.1 in Section 2, the complex Hessian comparison theorem on distance functions 
proved in [LW] and [CN] recently and the well-known heat kernel estimates of Li-Yau on 
complete Riemannian manifolds with nonnegative Ricci curvature, which states that 

C-\n) ( r 2 x (y)\ TTl . C(n) ( r 2 x (y) 
' exp — ^ < H(x, y, t) < exp ' 



V x (Vt) "V 3* ; " v V x (Vi) "V 5t 

for some C(n) > 0, where n is the real dimension of the manifold considered. As applica- 
tions we prove an elliptic 'monotonicity principle' for complex subvarieties in M. It can 
then be applied to prove a manifold version of Stoll's theorem. 

Before we prove a localized version of Theorem 2.1, we need to introduce some functions 
(notations). Let V be an analytic subvariety (of M) of complex dimension s. For the 
simplicity of the notation, Let H (x, y, r) be the fundamental solution to the heat equation 
— A) u(x,r) = 0. When x,y G V, we denote ( r Kr) rn ~ s H{x, y, r) by H v (x,y,r). For 
any fixed (x , t ) with x G V, we denote H v (x , y, t -t) by H( XOjto ) jV (y, t). For any p > 0, 
we also introduce a cut-off function 

V{x ,t ), P {Vi t )={ 1 ~2 J 

where f+(x) = max(/, 0) for any function /, r XQ (y) is the distance function (of M) from 
xq to y. It is easy to see that (p( Xo ,t ),p ^ s supported in B Xo (y/ p 1 — s{t — to)). 
The following simple lemma is useful. 

Lemma 4.1. On V, 

(4.2) (r A v)«.p(».^ - 

Here Ay denotes the Laplacian operator with respect to the induced Kahler metric on V 
(strictly speaking only regular part ofV). 

Proof. For any y G V, choose a complex coordinate [z\ • ■ • , z m ) as in the proof of Theorem 
2.1. Namely z a = on V for any a > s. We also use the index convention as in the proof 
of Theorem 2.1. Namely, 1 < i, j, k, • • • < s and s + 1 < a, b, c, ■ • ■ < m. Direct computation 
shows that 

>s- g lJ gc 3 
> 0. 

Here we have used the Hessian comparison theorem on the distance functions proved in 
[LW] (see also [CN, Corollary 1.1]). 

The following is a localized version of part (ii) of Theorem 2.1. 
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Proposition 4.1. Let 

(4-3) ^.so.to.tiC*) = / ¥ 3 (xo,ti),p(2/^)^(x ,to),v(2/^)^v- 

./V 

When in the right context we also it briefly denote by E\>. Then 
d f 

(4-4) j f E V (t) < -J IV-Llog^.t^.vlVtxo.tO.p^xo.toJ.V^V- 

Proof. The computation (2.5) implies that 

(4-5) ^ + A v ^j H {xo , toh v < -|V X log^ (xo , to)iV | 2 if (x0ito)jV . 

By Lemma 4.1 we have that 

dt Ev ^ = Iv {jft V{x ° M) 'p) H{*o,t ),V + <P(x ,ti),p (J^(*o,to),v) dA v 

+ ^( A V¥ 3 (xo,t 1 ),p)^(x ,t ),V - <P(x m),i> ( A vH( Xo , to ),v) dA v 

<- \V ± logH {x0ito)y \ 2 (p {xoM)jP H {x0jto)y dA v . 
Jv 



Here we have used the observation that 



/ , H ix0jto)y (V v <p {x0ttl)tP ,v)dS < 

Jd(vnB X0 (y/p2+s(t -t))) 



where v is the unit out- normal of d(V fl B Xo {^J p 2 + s(ti —t))) in V and dS is the area 
integral of d(V n B XQ ( ^ p 2 + s(h - t))). 

We denote the 2s-dimensional Hausdorff measure of set V fl B XQ (p) by Ay, Xo (p). As 
a consequence we have the following elliptic 'monotonicity principle' (can be viewed a 
Bishop-Lelong lemma on manifolds). 

Corollary 4.1. (Monotonicity principle) Let 8{s) = ^ 2 ^_ 4s . There exists C = C(m,s) 
such that for any p' G (0, S(s)p) 

(4 . 6) ^yif 1 <ci„i ^ M . 

Vx \P ) Vx \P) 
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Proof. The proof follows essentially the argument of Proposition 3.5 in [E2]. Applying 
Proposition 4.1 with t\ = t — S 2 (s)p 2 and to replaced by to + p' 2 . Using the heat kernel 
upper bound of Li-Yau we have that 

(4.7) 

C{m) 



H {xo!to+p/%v (y,t -5 2 (s)p 2 ) < (n(5 2 (s)p 2 + p' 2 )y 

(np 2 ) m - s 



V X0 (^/p' 2 + 5 2 (s)p 2 ) 



< C(m, s) 

< C(m, s) 



V X0 (6(s)p) 
{np 2 ) m - S 
V X0 (p) ' 



Here we have used the Bishop volume comparison in the last inequality. Notice that 
f(x ,ti),p( x ^° ~ d 2 ( s )p 2 ) — 1 an d supported inside B Xo (p). Hence (4.7) implies that 



(4.8) E v (t -S A (s)p')<C(m,s) f , 

Vx (P) 

By Proposition 4.1 we know that 

(4.9) E v (t - 5 2 (s)p 2 ) > E v (t - p' 2 ). 

On the other hand, by Li-Yau's heat kernel lower bound we also have that, for all y G 

B X0 (p% 

(4.10) 

Vx {P) 

Again we have used the Bishop volume comparison theorem. Notice that for y e B XQ (p') 

< + ' 2 W1 P' 2 + S(S 2 (S) P 2 + P' 2 ) 

(4.11) *o - p ) > i - — — 

> 1 -5 2 (s)(l + 2s) 
1 

~ 2" 

Combining (4.10), (4.11) we have that 

/2(m — s) a i ;\ 

(4.12) E v (t - p' 2 ) > C(m, s) P - V ( V ,\ Xo[P) . 

Combining (4.8), (4.9) and (4.12) we complete the proof. 
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Remark 4.1. In Proposition 3.1.1 of [Ml], a (considerably weaker) comparison on the 
relative volumes in the similar spirit of Corollary 4-1 was established for the zero divisors 
of holomorphic functions of polynomial growth, under further assumptions on M being of 
maximum volume growth and of quadratic curvature decay, using very different method. 

The following consequence of Corollary 4.1 is somewhat surprising. The result sharpen 
the Bishop-Gromov volume comparison theorem in the presence of compact subvarieties. 

Corollary 4.2. Let M m be a complete Kahler manifold with nonnegative holomorphic bi- 
sectional curvature. Suppose that M contains a compact subvariety V of complex dimension 
s. Then there exists C = C(m, s) > such that for S(s)p > p' >> 1, 

t 7" / \ / \ 2(m— s) 



<c ^ 



In particular, 

p^oo p2(m-s) 

An application of the above 'monotonicity principle' we can have another proof of The- 
orem 3.1 of [N4]. In fact the new proof gives a characterization of divisors defined by 
holomorphic functions of polynomial growth (also called 'polynomial functions' according 
to the notation in [W2]). 

Theorem 4.1. Let M be a complete Kahler manifold with nonnegative bisectional curva- 
ture. Let V be a analytic divisor of M . Define the Lelong number (elliptic) at infinity ofV 
by 

(A -\r>\ (v\ y 7r P 2 A^o(p) 

(4.13) ^oo(V) = sup hmsup — -. . . 

Assume further that H l (M, O*) = 0. Then V is defined by a polynomial function if and 
only i/^oo(V) < oo. Moreover, ifV = Z(f) for some f G P<i(M) (the space of holomorphic 
functions of polynomial growth with degree at most d) then there exists a C(m) such that 
for any xq G M , the Lelong number v(xq, V) at xq is bounded by 

(4.14) ^o,V)<C(m)i/oo(V) 
and 

(4.15) u x ( V) < C(m)d. 

Proof. Notice that (4.14) follows from Corollary 4.1 directly. First assume that V is the zero 
divisor of a polynomial function / (g Pd(M) for some d). We shall show that z^oo(V) < oo 
(in fact (4.15)). We follow the notation in Section 3 of [N4] (also Section 2 of this paper). 
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Let v(x,t) = J M H(x, y, t) A log \f\ 2 dv y and w(x,t) = -^v(x,t). By Corollary 2.1 we know- 
that 

tw(x,t) = (irt) / H{x,y,t) dAy. 
Jv 

By Li-Yau's lower bound estimate on heat kernel we have that 

Then (4.15) follows from (3.13) of [N4]. 

Now we assume that ^oo(V) < oo we prove that V is the divisor of a polynomial function. 
First, by the solution to Cousin problem II (directly from the vanishing of the co ho mo logy 
i? 1 (M, O*)) we know that there exists a holomorphic function / such that Z(f) = V. First 
we apply the 'moment type estimate" from [Nl] to estimate tw(x, t) from above. Note that 
w(x,t) = f v H(x,y,t) dAy(y) is well-defined due to the assumption that z^oo(V) < oo. By 
Corollary 2.1 we know that 

w{x,t)= I H{x,y,t)Alog\f\ 2 dfi{y). 

JM 

Applying the 'moment type estimate', Theorem 3.1 of [Nl] we have that 
(4.16) tw{x,t)<C{m)v OQ iy). 

Now we define v(x,t) = J* w(x, r) dr + log \f\ 2 (x). Then (4.16) implies that v(x,t) < 
C(m)^oo(V) log(t + 1) + log \ f\ 2 (x) for t ^> 1. One can also check that v(x,t) is a solution 
to the heat equation (J^ — A) v(x,t) = with v(x,0) = log |/| 2 (x). Now we apply the 
'moment type estimate' of [Nl] again we have that 

/ log \f\ 2 (y) dfi(y) < C(m)i/oo(V) log(r + 1). 

JB x (r) 

Now applying the mean-value inequality of Li-Schoen [LS] we conclude that / is of poly- 
nomial growth. 

Remark 4.2. In [St] (see also [Ru] for the codimension one case), Stoll proved that an 
analytic divisor V in M = C m is algebraic if and only if ^oo(V) < oo. The result was 
later generalized to the case M being an affine algebraic variety in [GK] by Griffiths and 
King. In fact, in [St] the result was proved for any analytic sets of C m . It is desirable to 
generalize our result to the high codimension case. 

The assumption on vanishing of the cohomology in the above theorem is satisfied, for 
example, when M is Stein (cf. Theorem 5.5.2 of [Ho],). 

In [N4], the author developed a parabolic approach to compare the vanishing order of 
a holomorphic function, at any fixed point, with the growth order at infinity. The method 
is sharp and effective. It turns out that the parabolic method there is also related to the 
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Nevanlinna theory for several complex variables. In order to illustrate this connection we 
need to recall some basic notations from the Nevanlinna theory (cf. [Gr]). 
We call that a function s : R + — > R+ has finite order if 

n m \ v logs(r) 

Urd(s) := hmsup < oo. 

r^oo logr 

For a / G 0{M), the space of holomorphic functions, we define the order of / in sense of 
Hadamard by 

Ord H (f) = Ord(log(A(r))) 

where A(r) = sup xeBo ^ \f\(x) with o G M being a fixed point. Following [Gr, St] for any 
analytic subvariety of complex dimension s we define 

nvM= V„(r) 

and 

f r dr 
N v (x ,r)= / (n v (x ,r) - n v (x o ,0)) h n v (x , 0) logr. 

A^v(^o 5 ?") ( n v(xo, r)) is called the counting function since for m = 1 and V being the zeros 
of a holomorphic function, it simply counts the number of zeros in B(xq, r). If Xq does not 
lie in V one has Ny(xo, r) = J Q r ny(xo, t)^. One can view the estimate (4.15) as bounding 
the counting function by the growth order, a Nevanlinna type inequality. The following 
result is a further generalization of (4.15), whose correspondence in the Euclidean space is 
also known as the transcendental Bezout estimate. 

Corollary 4.3. Let M be a complete noncompact Kdhler manifold with non-negative Ricci 
curvature. Let f G O(M) be a holomorphic function of finite order. Let Z(f) be the zero 
divisor of f . Then 

(4.17) Ord(N z (x , r)) < Ord H (f). 



Proof. The proof follows the same line of argument as the proof of Theorem 4.1. We leave 
it to the interested readers. 

We found the connection between the parabolic equations (as well as the related differ- 
ential Harnack inequality) and the Nevanlinna theory quite interesting. 

§5 Heat kernel and the reduced volumes. 

In this section we consider the fundamental solution to the time-dependent heat equation 
with metrics deformed by Kahler-Ricci/Ricci flow (1.1). As we explained in Section 2 that 
the fundamental solution H(x,y,t,to) satisfies (1.2). We shall prove a sharp lower bound 
on H(x,y,t,to). In [N5] we derived an estimate for the fundamental solution satisfying 
the time-dependent heat equation itself (instead of (1.2)). The estimate there is only valid 
for short time interval. In [Gu], a rough lower bound was obtained through the earlier 
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mentioned Harnack inequality, Theorem 1.3. Notice that the result in [Gu] is not sharp 
and the result in [N5] is only sharp in the exponents. In this section we show a sharp lower 
bound on the heat kernel in the case that M is either a complete Riemannian bounded 
nonnegative curvature operator or a complete Kahler manifold with bounded nonnegative 
bisectional curvature. For the sake of simplicity we only state the result for the Kahler- 
Ricci flow and leave the Riemannian analogue to the interested readers. Before we state 
our result we need to recall some notations and computations from [FIN] (which follows 
closely the computation in [P]). For the simplicity we assume that to = 0. 

Let g(t) solve the Kahler-Ricci flow on M m x [0, T] (m = dimc(M) and n = 2m). Fix 
xq and let 7 be a path (x(r]),r]) joining (xo, 0) to (y,t). Following [P] and [LY, FIN] we 
define 

(5.1) £+(7)= f Vv(K + 4h'(v)\ 2 ) d V . 

Jo 

Let X = "f'(t) = dz d ^ of^- and let Y be a variational vector field along 7. Here |7'(^)| 2 = 

9a8 dZ at^ ^df* ■ Using C + as energy we can define the £ + -geodesics and we denote L + {y, t) 
to be the length of a shortest geodesies jointing (xo,0) to (y,t). We also define 

Following the first and second variation calculation of [P] (see also [FIN]) we compute that 



(5.2) | V £ + | 2 = -^+^+ + A 



t t 3 / 2 ' 



(5.3) «t= W - K 



' + 



dt 2t 3 / 2 t 

n K 

(5.4) A£+ < K + ^ 



2t 2t 3 / 2 ' 



Here 



K := I r] 3/2 H(X)dm 
Jo 



where H{X) := 811/ dt + 2(VTl, X) + 2(X, VIZ) + 4Ric(X, X) + TZ/t is exactly the traced 
Li-Yau-Hamilton differential Harnack expression in [CI] applying to the (1,0) vector field 
2.V. 

Theorem 5.1. Let (M m , g(t)) be a complete solution to Kahler-Ricci (Ricci) flow with 
bounded nonnegative bisectional curvature (curvature operator). Let H(xo,y,t) be the fun- 
damental solution to the time dependent heat equation (satisfying (1.2) and H(xo,y,0) = 
S Xo (y)). Then 

u{x, t) := j—^ exp (-£+(x, t)) 
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(5.5) ( ^- - A -11 ) u(x,t) < 0. 



satisfies 

! .< :/ 

In particular, 

(5.6) < H(x ,y,t) 

and 



e^°'°\t):= [ u{x,t)dm 

J M 



is monotone decreasing. Moreover, the equality in (5.5), or (5.6) implies that M is a 
gradient expanding soliton. 

Proof. First (5.2)-(5.4) implies that 

| - A - *) ((J,;; exp (-M., *))) = "f exp (-* +to , *))) < 0. 

Here we have used fact that K > under the assumption that M has bounded non- 
negative bisectional curvature. Also if the equality holds it implies that K = 0. This 
further implies that M is an expanding soliton from the computation in [FIN] . In order to 
prove (5.6) one just need to apply the maximum principle (cf. [NT2, N5]) and notice that 
lim^o ^\) m ex P(~^+(^)) = 3x (y)- The equality case follows from the consideration of 
the equality in [FIN]. 

Remark 5.1. In [CY], Cheeger and Yau proved that for complete manifolds with nonneg- 
ative Ricci curvature, the heat kernel H(x, y, t) has the lower bound estimate: 

H(d(x,y),t) := (^exp(-^|^) < H(x,y,t) 

by showing that the transplant of the Euclidean heat kernel H(d(x,y),t) is a sub-solution 
to the heat equation. Recall that in [P], Perelman first discovered that there is a similar 
result for the backward Ricci flow, even without any curvature sign assumptions ( this is the 
astonishing part of Perelman's work). Namely he proved, for -§^g%j = 2Rij and the reduced 
distance £(y,r), defined formally by the same expression as £ + , that 

(5.7) (J^-A + nju(y,r)<0 

where ^ 

u{y,r) = — — -^exp(-£(y,r)). 

(47TT) 2 

In particular, Perelman further showed (in Corollary 9.5 of [P]) that u gives a lower 
bound for the heat kernel of the backward heat equation ((J^ — A)v = 0). (Notice 
that the heat kernel H(x, y, r) for the backward heat equation satisfies the Schrddinger 
eqaution (J^ — A + H(x, y, r) = 0. One should also refer to [LY, Theorem 4.3] for 
a precedence of Perelman's second variation computation in [P, Section 7].) Thus both 
Perelman's monotonicity of the reduce volume and Theorem 5.1 above can be viewed as 
nonlinear analogue of the earlier work of Cheeger- Yau in [CY] and Li- Yau in [LY] . 

Tracing the computation of [FIN] we also have the following estimates for u{x,t). 
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Proposition 5.1. Assume that (M,g(t)) be a Kahler-Ricci flow with bounded nonnegative 
bisectional curvature on M x [0, T). Let u(x,t) be as above. Then 

(5.8) log(w) Q ^ + R a p + -g a p > 0. 

The equality holds if and only if (M, g(t)) is an expanding Kahler-Ricci solition. 

Notice that (5.7) is equivalent to the estimate (1.3) on u(x, t) a positive solution to the 
forward conjugate heat equation. On the other hand u(x,t) is only a sub-solution to the 
forward conjugate heat equation. 

Remark 5.2. In [FIN, Corollary 2.1] we showed that under the assumption of bounded 
nonnegative curvature operator there exists a same estimate as (5.8) (without the bars). 
This suggests that one may have the similar estimate as Theorem 1.2 for the Riemannian 
case on positive solution u of the forward conjugate heat equation. Indeed, this can be 
shown if one proves, under the assumption of M having bounded nonnegative curvature 
operator, that 

1 

(5.9) Yij := V i V j K-2VkRijV k {logu) + 2R ikjl Vk(logu)Vi{logu) + R ik R jk + -R lJ > 0. 

This claim is based on the following computation which is essentially due to Chow. Note 
that = on gradient expanding solitons and its trace is the same as the trace of 
Hamilton's matrix Li- Yau- Hamilton expression in [HI]. 

Lemma 5.1. Let 

Zij = R l3 + {logu)ij + — gij . 

Then 
(5.10) 

^ - A L ^j = Y^ + 2V fc Z ij -V fc (logu) 

+ Z ik (\og(u) jk - R jk - ^.9jkj + (log(u) ik - R ik - Y t Sl ^j 

where AlT}^ = Arjij + 2Ri k j\r\ k \ — RikVkj — RikVik ^ s ^he Lichnerowicz operator acting on 
the symmetric 2-tensor rjij. 

The above computation was first carried out in [Ch2] for the backward Ricci flow, with 
Y^ is replaced by Hamilton's matrix Harnack expression for shrinkers. Notice also that 
the matrix Harnack expression Y^ and Hamilton's expression for Ricci expanders are the 
same if the manifold is Kahler. 

One can obtain some upper bound on the heat kernel H(x,y,t,0) using the Harnack 
inequality proved in Corollary 2.1 and the fact that J M h(x, y, t, 0) dfit = 1. But the result 
is not as satisfactory as for the fixed metric case of Li- Yau. Hence we shall leave this 
to a later investigation. On the other hand, the localization technique of Ecker can be 
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applied to the reduced volume for the Ricci expanders (defined in [FIN]) to obtain the 
monotonicity of a localized reduced volume without assuming bisectional curvature (nor 
curvature operator) being nonnegative. Recall that in [FIN], the authors proved that, if M 
is a closed manifold, the forward reduced volume 

e { *°'°\t) := / u(x,t)dnt 

J M 

is monotone non-increasing, where u(x,t) = e ^ v ^ m ) . Here we define u for the Kahler case 
to be coherent with our previous discussions. This monotonicity has severe restriction since 
the reduced volume 9+'°\t) is only meaningful when M is compact. We shall show the 
monotonicity of a localized reduced volume, which is well-defined on complete manifolds, 
to compensate such restriction. In order to put the result in its general form we denote 
by /^ c °' t °- ) the forward reduced distance with respect to (xo, to) if one replaces the reference 
point (x ,0) by (xo,to). Correspondingly we denote by u^ ^ anc [ q(*°> *o) j the reduced 
volume function and the reduced volume with respect to (xo,to), respectively. Let 



and 



Then (5.2)-(5.4) imply that 

(5.11) (J^_+A-Kj u( Xo > to \x,t) <0 



and 
(5.12) 



d_ 
dt 



-A 



(xiM) 



(x,t) < 0. 



Notice that if IZ is uniformly bounded from below, (ft^p 1 ^ {x,t) is compactly supported. 
We then have the following localized monotonicity of the reduced volume. 

Proposition 5.2. 

(5.13) ±e^\t)<0 
where 

(5.14) e^ t0 \t)= [ ^\x,t)u^\x,t)d^ t . 

J M 
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It turns out that one can also construct a subsolution to the heat equation with compact 
support, similar to (p^^^x^t) defined above. This is through Perelman's reduced distance 
^(*i,*i)( a;j r ) ? w ith r = t 1 -t. (^'^(Mi - 1) is only defined for t < t x . Please refer to 
[P, Section 7] for the detailed discussions.) Here we use (xi,t±) to specify the reference 
space-time point with respect to which the reduced distance is defined. Recall that from 
[P, Section 7], L^ Xl,tl \x,t) := (t\ — t)v- Xl,tl \x,t\ — t) (again we do not have the factor 4 
since we are in the Kahler setting) satisfies the differential inequality: 

Then in the case 1Z is bounded from below, we may define a compact supported function 

CM^(i- p " WM ; B "- tl ) ; 

It is easy to see that 

(5-15) (l-^^\ x ,t)<0. 

This gives another localization on the forwarded reduced volume. 
Proposition 5.3. 

(5.16) j/^/°\t)<Q 
where 

(5.17) efg*>\t)= [ iP^\x,t)u^{x,t)diH. 

J M 



The similar idea can also be applied to Perelman's entropy. Let (xo,to) ^ e a fi xe d 
space-time point. Let t = to — t and u(x,r) be the fundamental solution of the backward 
conjugate heat equation — A + 1Z centered at (xo,to). Write u = ^ T ^ m and define 
v = (t (2Af - |V/| 2 + TZ) + f - 2m) u. Then it was proved in [P] that 

(5-18) (J^-A + k)v=-t (\R aB + V a V- p f - ^g a p\ 2 + \V a V p f\^J u 

and v < 0. Applying (5.12) or (5.15) we can have the following local monotonicity formulae. 
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Proposition 5.4. 

(5.19) 

j t [J <- J T (\ R c*+ v - v ^/ - ^m 2 + i v « v /3/i 2 ) ^b ti} ^ 

and 
(5.20) 

j t (/ m ^) ^~j M T (l^a/S + V a V^/ - 1^| 2 + |V a V^/| 2 ) 

Proof. The proof is just direct computations and integration by parts, using (5.12), (5.15) 
and (5.18). 

There is no difference for the Riemannian cases. The formulae are exactly the same 
except some factors caused by the definition of the operators. 

We conclude this section by an application of Theorem 5.1 to the study of the large 
time behavior of Kahler-Ricci flow. Let (M m , go) be a complete Kahler manifold with 
bounded nonnegative holomorphic bisectional curvature. If we further assume that (M, go) 
has maximum volume growth, namely for any x, the volume of the ball B(x,r), V x (r) is 
bounded from below by Sr 2m for some positive 5, in [N6] we proved that the Kahler-Ricci 
flow (1.1) has long time solution with the initial data g(x,0) = go(x). Moreover, there 
exists a constant A = A(M) > such that t ^> 1, 

(5.21) tTl{x,t)<A. 

Applying Theorem 5.1 to this situation we can show the following result. 

Corollary 5.1. Let (M,go) be a complete Kahler manifold as above. Let g(x,t) be a 
long time solution to Kahler-Ricci flow as above. For any (xj,tj) with tj — > oo, define 
gj(t) = j-g(tjt). Then the pointed sequence (M,Xj,gj(x,t)) sub-sequentially converges to 
a gradient expanding Kahler-Ricci soliton (M^, x oo,<7oo(£))- 

Proof. For any t, denote by B t (x,r) the ball of radius r with respect to g(t), and by 
Vt(x, r) the volume of this ball (with respect to g(t)). By Theorem 2.2 of [NT3], we know 
that Vt(x,r) > 5r 2m . Together with (5.21) we have the injectivity radius bounded from 
below uniformly for gj(t). Therefore, by Hamilton's compactness theorem, (M 7 xj, gj(t)) 
subsequentially converges to (M^, x^, goo(t)), a solution to Kahler-Ricci flow defined on 
Moo x (0, oo). The only thing we need to prove is that (Moo, 9<x>(t)) is an expanding soliton. 
It is easy to see that (Moo,<7oo) also has bounded nonnegative bisectional curvature and 
the maximum volume growth. By Corollary 1 of [N6], we know that M^ is topologically 
R 2m . The fact that it is an expanding soliton follows from the monotonicity consideration. 
Let us adapt the notations from Theorem 5.1. Let (x ,0) be a fixed reference point, with 
respect to which we have the forward reduced distance function £+(x,t), and the second 
forward reduced volume function u(x, t) and the second forward reduced volume 6+(t). By 
the proof of Theorem 5.1 we have that 
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from which we can conclude that 

r>2t. 



Hence 



§+(21,) - h(tj) > f ' KudfH dt. 



Taking limit we have that for the limit metric 

r-2 



1 



if« (it = 0. 



The conclusion now follows from [C2] or Theorem 4.1 of [N2] since K = implies that the 
linear trace LYH quantity H achieves its minimum (zero) somewhere. 

Remark 5.3. Corollary 5.1 simply says that the blow-down limit of a Type III solution 
is an expanding soliton. The result is in some sense dual to Proposition 11.2 of [P]. It is 
also similar to the situation studied in [Hu] for the mean curvature flow. The maximum 
volume growth assumption can be weaken to a certain K-noncollapsing condition in terms 
of the lower bound of 9 + (t) . 

It has been proved in [CT1] (see also [CT2] for further more recent progresses) that a 
gradient expanding Kahler soliton must be biholomorphic to C m . 

One can have a similar result for the Riemannian case, replacing the nonnegativity of 
the bisectional curvature by the nonnegativity of the curvature operator. 

§Appendix: A parabolic relative volume comparison theorem. 

Recall that Cheeger and Yau proved that on a complete Riemannian manifold with 
nonnegative Ricci curvature, the heat kernel H(x, y, r) (the fundamental solution of the 
operator (J^ — A)) has the lower estimate 



(Al) H(x,y,r) > — — -^exp 



r 2 (x,y) 



(Aut) 2 \ At 

where r(x,y) is the distant function on the manifold. This fact can be derived out of the 
maximum principle and the differential inequality 

Integrating on the manifold M, this differential inequality also implies the monotonicity 
(monotone non-increasing) of the integral 

(A2) V M :=[ expf-^ 
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In [P] , Perelman discovered a striking analogue of this comparison result for the Ricci flow 
geometry. More precisely, he introduced a length function, called the reduced distance, 

4 ,s(t) (y, f) := inf J Vrh'(r) | 2 dr 

for all 7(r) with 7(0) = xo, 7(t) = y (in the right context, we often omit the subscript 
xo,g(r)), and a functional, called the reduced volume, 

(A2') V 9{t) (x ,t):= [ — i-^exp(-%,r)) 

with respect to a solution gij(x, r) to the backward Ricci flow 

7^: #17 = 2-Rij 

on M x [0,a]. Perelman proved further that V g ^ T )(xo,r) is monotone non-increasing in r. 
In fact, he showed this monotonicity result via a space-time relative comparison result and 
used it to give a more flexible proof of the K-noncollapsing on the solution to Ricci flow, 
defined on M x [0, T], for any given finite time T with given initial data. 

In [P] , Perelman also discovered a more rigid monotone quantity, the entropy functional. 
These two functionals can be related through a differential inequality of Li-Yau-Hamilton 
type. One can refer to [CLN] for an exposition on this relation. The entropy is stronger but 
less flexible. The entropy monotonicity has its analogue for the positive solutions to linear 
heat equation on a fixed Riemannian manifold with non-negative Ricci curvature. This 
was derived in [N3] , where the author also observed that the entropy defined on a complete 
Riemannian manifold M with non-negative Ricci curvature is related to the volume growth 
of the manifold. More precisely, let 



(A3) W(/,t):=/ (T|V/| 2 + /-n) 7 -^ 3r dp 

for any r > and C 1 function f with f , 6 1 n da = 1. If u(x, r) = e 1 n is a solution 

to the heat equation, it was proved in [N3] that W(/, r) is monotone non-increasing in r. 
Moreover, it was shown that 

lim W(/,r)=log(lim V Xo (r)) 
if u(x, r) is a fundamental solution originated at x$. In fact, 

lim V X0 (r)= lim 
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where V Xo (r) is the volume of ball of radius r, which is independent of xq. (The limit on 
the right hand side above is also called the cone angle at infinity). 

Motivated by this close connection between the Ricci flow geometry for a family of met- 
rics and the Riemannian geometry of a fixed Riemannian metric, it is nature to seek a 
localized version of the above mentioned Cheeger-Yau's result on the monotonicity of the 
reduced volume V Xo (t). In fact, a local version of the heat kernel comparison was also 
carried out in the paper of [CY] from a PDE point of view. However, Perelman's localiza- 
tion in the case of Ricci flow geometry is along the line of comparison geometry, and very 
much different from Cheeger-Yau's localization consideration (by considering the Dirichlet 
or Neumann boundary value problem for the heat equation). This leads us to formulate 
a new relative (local) volume comparison theorem in this section, which is the linear ana- 
logue of Perelman's formulation for Ricci flow. It can be viewed a parabolic version of the 
(classical) relative volume comparison theorem in the standard Riemannian geometry. The 
interested readers may want to compare these two results in more details. The relation 
between the classical relative volume comparison theorem and the new one here is very 
similar to the one between the monotonicity formula for the minimal submanifolds in M. n 
and Huisken's monotonicity for mean curvature flow in M n (as well as Ecker's localized 
version) . 

In order to state our result, let us first fix some notations. It is our hope that the 
exposition below is detailed enough to be helpful in understanding [P] better. Fix a point 
xo G M. Let 7(7-) (0 < t < f) be a curve parameterized by the time variable r with 
7(0) = xq. Here we image that we have a time function r, with which some parabolic 
equation is associated. Define the £-length by 



Jo 

We can define the £-geodesic to be the curve which is the critical point of £(7). The simple 
computation shows that the first variation of C is given by 



where Y is the variational vector field, from which one can write down the £-geodesic 
equation. It is an easy matter to see that 7 is a £-geodesic if and only if 7(a) with a = 2y / r 
is a geodesic. In another word, a £-geodesic is a geodesic after certain re-parametrization. 
Here we insist all curves are parameterized by the 'time'-variable r. One can check that 
for any v E T Xo M there exists a £-geodesic 7(7") such (7(0")) |<r=o = v - Notice that 
the variable a scales in the same manner as the distance function on M. So it is more 
convenient to work with a. 

We then define the ^-exponential map by 



(A4) 




(A5) 



5£( 1 ) = 2Vf(Y,X)(f)-2 f T ^((V X X+^X,Y)) dr, 



Cexp v (a) := 7„(tr) 



if 7^(0-) is a £-geodesic satisfying that 



Bm-( 7 ») = U . 
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It is also illuminating to go one dimensional higher by considering the manifold M = 
M x [0, 2y/T] and the space-time exponential map exp(5°) = (£exp<^ (a), a), where v a = 

a 

(v a ,a). Denote ^- simply by v 1 and (V,l) by v 1 . Also let 750(77) = exp(r]v a ). It is easy 
to see that 

This shows that 

<iexp|( 0i o) = identity. 
Computing the second variation of £(7) gives that 

(A6) 5 2 £( 7 ) = 2y/T(V Y Y,X)+ [ 2^{\V X Y\ 2 -R(X,Y,X,Y)) dr 

Jo 

where Y is a given variational vector field. Let (y,(r) = (£exp v (a),a). Consider the 
variation which is generated by (w, 0). Namely consider the family 

U(s,a) = (Cexp v+sw (a),a) = exp((a(v + sw), a)). 

Direct calculation shows that 

^(o,o) = Ki), ^(0,0) = 

and that the Jacobi field J w (a) = (J w (a),0) is given by 

-fe\*=o( a ) = dexp((crw,Q)) 

with the initial velocity 

X7 « I 

ds 



Va l^| 4=O )(0)=K0). 



(One can define the Jacobi operator to be the linear second order operator associated with 
the quadratic form in the right-hand side of (A6). One call a vector field along 7 a C- 
Jacobi field if it satisfies the Jacobi equation. It is easy to show that the variational vector 
field of a family of £-geodesics satisfies the Jacobi equation as in the standard Riemannian 
geometry. In fact, the C- Jacobi field turns out to be just the regular Jacobi-field after 
re-par ametrization.) This shows that 

(A7) (d£exp) w (a) = J w (a) 

and 

(A8) dexp((aw, 0)) = J w (a). 

Now we can conclude that (y, a) is a regular value of the map exp (and y is a regular value 
of £exp(.)(cr)) if and only if that any Jacobi field J with initial condition as above does 
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not vanish at a. We can introduce the concept of conjugate point (with respect to xq) 
similarly as in the classical case. We can define the set 

D(a) C T Xi) M 

to be the collection of vectors v such that (£exp v (<x), a) is a ^-geodesic along which there 
is no conjugate point up to a. Similarly we can define the set 

C(a)CT X0 M 

to be the collection of vectors v such that (C v (a),a) is a minimizing £-geodesic up to 
a. One can see easily that D(cr) and C(a) decreases (as sets) as a increases. For any 
measurable subset A C T XQ M we can define 

D A {a) = A(~)D(a) and C A {a) = A n C{a). 

Following [P] we also introduce the ^-'distance' function. 

Zx (y,f) = —=L XQ {y,f), where L Xo (y,f) = inf £(7)). 

Here our I is defined for a fixed background metric. We also omit the subscript xq in the 
context where the meaning is clear. The very same consideration as in [P], as well as in 
the standard Riemannian geometry, shows that 

(A9) |W| 2 = -£ 

T 

(A10) t T = --£ 

T 

and 

n 1 r 3 

(All) A£<—-— T*Ric(X,X)dT 

2r 7-2 Jo 

where X = Y(t) with 7(r), < r < f being the minimizing £-geodesic joining xq to y. 
Putting (A9)-(A11) together, one obtains a new proof of the result of Cheeger-Yau, which 
asserts that if M has nonnegative Ricci curvature, then 

— l(y .t) 

Namely — is a sub-solution of the heat equation. In particular, 



(4-kt)~ 



dr J M (4ttt)2 
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Using the above geometric consideration, one can think the above result of Cheeger-Yau as 

— £(y,r) 

a parabolic volume comparison with the respect to the positive measure w ■ dfi. Recall 

(47TT)"? 

that the well-known Bishop-Gromov volume comparison states that if M has nonnegative 
Ricci curvature 



d / 1 

dr \ r n 



dA | < 

S* (r) 



where S XQ (r) denotes the boundary of the geodesic ball centered at xq with radius r, 
dA is the induced area measure. The by-now standard relative volume comparison can 
be formulated in the similar way as above. Let A be a measurable subset of >S n_1 C 
T Xo M one can define CU(r) to be the collection of vectors rv with v G A such that the 
geodesic exp XQ (sv) is minimizing for s < r, where exp (•) is the (classical) exponential 
map. Then the classical relative volume comparison theorem (cf. [Gr]) asserts that if M 
has nonnegative Ricci curvature 




dA | < 0. 

exp{CA(r)) 



The following is a parabolic version of such relative volume comparison theorem parallel 
to Perelman's work on Ricci flow geometry. 

Theorem Al. Assume that M has nonnegative Ricci curvature. Then 
(A13) — / - dfjL < 0. 

Proof. The proof follows the similar argument as in [P]. Using the notation in the above 
discussion, first we observe that 



r e -t(y,r) r e -t(v,r) 

/rex PCAM (T) (47rr)f ^ J Ca ( t ) (47rr)f 



^J(t) d/io, 



where J(r) is the Jacobian of £exp(.)(r). Since Ca(t) is decreasing in r, it suffices to show 
that 

d / e -e(y,r) \ 
< A14 > Tr {T^W J(T] ) £ °' 

This follows from (A9)-(A10), the fact that 

A e -'=(-^ + W X)(r)) e -« = 

(since V£(f) = A), and the claim that 

d n 1 r f s 

(A15) f-logJ(f) < - - / riRic(A,A)rfr. 

«t 2r r Jo 
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The estimate (A15) follows exactly the same proof as in [P], via the second variation 
formula (A6) and its consequence on the Hessian comparison: 

Hess(L)(y,y) < f 2V^{\V X Y\ 2 -R(X,Y,X,Y)) dr 
Jo 

where Y(r) is a vector field along the minimizing £-geodesic 7(7") joining xq to y with 
7(f) = y, satisfying Y(f) = Y and 

V X Y-—Y = 0. 



Better comparison between Theorem Al and the classical relative volume comparison 
can be seen by comparing M with M, M x {a} with S Xo (a). Notice that in [P], one 
does need such a result (not just the global version on the monotonicity of V g ^(xo, r)) to 
prove the non-collapsing result on finite time solution to Ricci flow. This and the potential 
application to the study of Riemannian geometry of our linear version, justify the spelling 
out of this result in spite of its simplicity. 

Remark. One can formulate the similar parabolic version of the volume comparison result 
for the case Ric(M) > —{n — 1)K . We leave that to the interested readers. 

In [N3], we showed that the lower bound on the entropy info< T <T M r ) (please see [N3] 
for definition) implies the non-collapsing of volume of ball of radius r for r 2 < T. This 
corresponds (but is considerably easier than) the K-non-collapsing result of Perelman on 
the solutions to Ricci flow, via the monotonicity of entropy functional. Using Theorem Al, 
one can have the following consequence, which is just the linear version of the second proof 
of Perelman on his celebrated K-non-collapsing result. 

Lemma Al. Let M be a complete Riemannian manifold with nonnegative Ricci curvature. 
Assume that B(xq, r) is K-collapsed in the sense that 

(A16) V(x ,r)<Kr n . 

Then there exists C = C(n) > 

(A17) V(K^r 2 ) < C ( v^ + exp ( ^ 

Namely, the smallness of the relative volume V ^°' r ' ) is equivalent to the smallness of the 
'reduced volume' V{r 2 ). 

Notice that the lemma can also be proved by direct computations and the Bishop volume 
comparison theorem, without using Theorem Al. 
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